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Abstract. In this paper, we revise Maxwell's constitutive relation and formulate a 
system of first-order partial differential equations with two parameters for compressible 
viscoelastic fluid flows. The system is shown to possess a nice conservation-dissipation 
(relaxation) structure and therefore is symmetrizable hyperbolic. Moreover, for smooth 
flows we rigorously verify that the revised Maxwell's constitutive relations are compatible 
with Newton's law of viscosity. 



1. Introduction 

Maxwell fluids are among macromolecular or polymeric fluids. A large number of 
experiments indicate that polymeric fluids exhibit elastic as well as viscous properties [T|. 
Thus, they are quite different from small molecular fluids. The latter have viscosity as 
the main feature, are satisfactorily characterized by Newton's law of viscosity 

Vv + {Vv) T - - V ■ vl] - reV ■ vl, (1.1) 
3 

and are also called Newtonian fluids. Here r = r(x, t) is the stress tensor of the fluid 
at space-time (x,t), v is the shear viscosity, re is the bulk viscosity, v = v(x,t) is the 
velocity, V is the gradient operator with respect to the space variable x = (xi, X2, X3), the 
superscript T stands for the transpose operator, and I denotes the unit matrix of order 
3. Combining Newton's law of viscosity with the conservation laws of mass, momentum 
and energy, one gets the classical Navier-Stokes equations. 

To account for the elastic properties of polymeric fluids, Maxwell combined Newton's 
law of viscosity with Hooke's law of elasticity and proposed the following constitutive 
relation [8] 

2 

€Tt + r = -u\Vv + (Vv) T - -V • vl] - reV • vl. (1.2) 

3 

Here e is the ratio of the viscosity and the elastic modulus. A Maxwell fluid is that obeying 
the constitutive relation (11. 2p . This relation reflects that the stress tensor responds to 
the fluid motion in a delayed, instead of instant, fashion. It has motivated many more 
realistic and nonlinear constitutive relations, including the well-known upper-convected 
Maxwell (UCM) and Oldroyd-B models [3]. 

In this paper, we revise Maxwell's constitutive relation (II. 2p . combine the conservation 
laws and formulate the following partial differential equations for compressible viscoelastic 
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fluid flows: 

d tP + v-{pv) =0, 

dt(pv) + V • {pv ® v + pl) + f V • ri + f Vr 2 = 0, 

(1.3) 

<9 f ri + - [Vu + (Vv) T - | V • i>l] = 



a(T 2 + iv-D = -^v 

Here p is the density of the fluid, ® denotes the tensorial product, p = is the 
hydrostatic pressure, t\ is a tensor of order two, e\ and e 2 are two positive parameters, 
and t 2 is a scalar. This is a system of first-order partial differential equations, with domain 

G := {(p,pv,T U T 2 ) : p > 0}. 

In ( II. 3p there are 14 equations (for three-dimensional problems). Note that the [•••]- 
term in the Ti-equation is symmetric and traceless. It is easy to see that n is symmetric 
and traceless if it is so initially. When n is symmetric and traceless, the number of 
independent equations in (JOJ) reduces to n — 10. Throughout this paper, we assume 
that Ti is symmetric and traceless. 

We will show that the first-order system (JOJ satisfies the entropy dissipation condition 
proposed in [13J. This particularly implies that the system is symmetrizable hyperbolic. 
Moreover, we will show that the revised Maxwell's constitutive relations (the Tp, t 2 - 
equations) in (II. 3p are compatible with Newton's law of viscosity (II. ip for small ei and 
e 2 . To see this, we rewrite the two r-equations in (JOJ as follows and iterate them once 
to obtain 

ti = -eiu[Vv + {Vv) T - |V ■ vl] - ejvdtn, 

= -e 1 u[Vv + (Vv) T - |V-vJ] +0(4), 
r 2 = -e 2 KV • v - elnd t T 2 , 

= -e 2 nV -v + 0(el). 

Substituting the truncations into the momentum equation in (II. 3p . we obtain the classical 
isentropic Navier-Stokes equations. In this sense, Newton's law (II. ip is recovered. 

A major part of this paper is devoted to a rigorous justification of the compatibility 
above. To do this, we employ the convergence-stability principle [T2| [2] for initial-value 
problems of symmetrizable hyperbolic systems and prove that, as e\ and e 2 go to zero, 
smooth solutions to the first-order system exist in the time interval where the isentropic 
Navier-Stokes equations have smooth solutions and converge to the latter. Namely, we 
show that the first-order system (11.31) is a diffusive relaxation approximation to the isen- 
tropic Navier-Stokes equations. 

Let us remark that, despite being quite similar, the present problem is very different 
from those studied in J5j [H]. In fact, when writing (JOJ) in its quasilinear form, the 
coefficients of - in the left-hand side depend on p and pv . Therefore, our problem does 
not possess the parabolic structure required in [S]. It also differs from that in [TT] . for 
p and pv are not dissipative quantities. Because of these, our analysis contains some 
innovative treatments relying on the specific structure of (JOJ). 
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We end this introduction by mentioning some other related works known to the author. 
Studying the motion of complex fluids involves many challenging and interesting partial 
differential equations [10] and has attracted much attention in recent years (see [7] and 
references cited therein). Most mathematical literature are concerned with well-posedness 
of incompressible flows governed by partial differential equations with upper-convected 
derivatives included. It seems that there are very few results on the Newtonian limit of 
non- Newtonian fluid flows. The only one known to this author is [S] , which was concerned 
with incompressible viscoelastic fluid flows of Oldroyd type for solutions in the Besov 
spaces. 

The paper is organized as follows. In Section 2 we show that the first-order system 
(11. 3p satisfies the entropy dissipation condition proposed in [13]. Section 3 is devoted to 
a precise statement of our compatibility result. A key error estimate is derived in Section 
4. 

2. Entropy dissipation structure 

In this section, we show that the first-order system (11.31) satisfies the entropy dissipation 
condition proposed in [13]. To this purpose, we define 



z z 



and compute, for smooth solutions to ( 11 .3D . 

d t (^(p) + 2p\v\ 2 + 2rl + \r 1 \ 2 ) 



(2.11 



+ V ■ (4$(p)u + 2p\v\ 2 v + 4pv + 4^ + 4^) = -g| - 

Here |ti| 2 is the trace of the matrix t^T\ and we have used the fact that T\ is a symmetric 
and traceless tensor. It is easy to verify that 

v = V (U) := 4$(p) + 2pM 2 + 2r 2 2 + |tx| 2 

is a strictly convex function of U := (p, pv , r 1; r 2 ) T , provided that the pressure p = p(p) 
is strictly increasing with respect to p > 0. Thus, the first-order system (11. 3ft fulfills 
the entropy dissipation condition in [13]. This particularly implies that the system is 
symmetrizable hyperbolic. 
Set 

w = (p,pv) T and z = (t 1 ,t 2 ) t . 
We may rewrite the first-order system (I1.3P (with e\ = e 2 = e for simplicity) as 

d t w + V ■ fj(w) X] + \ CjZ x = 0, 

(2.2) 

9tz + U2j9j{w) Xj =-±Sw, 

where Cj is a constant matrix and S = diag(^ _1 /9, k^ 1 ). Moreover, from the form of 
j] = i](U) and the strict convexity it follows that rj ww (w, z) is a symmetric positive-definite 
matrix, rj zz (w,z) is a constant diagonal (and positive-definite) matrix, r] ww (w, z)fj W (w) is 
symmetric, and 

r} ww (w, z)Cj = g jw ( w ) Tr lzz(w, z) T . (2.3) 
The last two statements are based on (12. ip . 
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3. Compatibility theorem 

This section is devoted to a precise statement of our compatibility result. For the sake 
of simplicity, we assume that ti = e 2 = e in what follows. 

Let p = p(x, t) and v = v(x, t) be the density and velocity of the Newtonian fluid. Then 
they obey the conservation laws of mass and momentum 

d t p + v-{pv) =o, 

d t {pv) + V • {pv <8> v + p{p)I) + V ■ r =0 
together with Newton's law of viscosity ( 11. li 



Vv + (Vv) T - -V ■ vl 



— • vl. 



Namely, they solve the isentropic Navier-Stokes equations. 
Define 

2. 

p e = p, v e = v, r lt = -eu 



Vv + (Vv) 1 - -V • vl 



T 2 e 



-e«;V • v. 



We have 



d t pe + V • (p e v e ) 

d t { Pe v e ) + V • {p e v e ®v e + p{p e )I) + iV • t u + 7Vr 2e 



0, 
0. 



d t r u + i + (Vv e 



i T - §V • u e J 



. t 2e 
Ke 2 



(3.1) 



Our compatibility result can be stated as 



Theorem 3.1. Suppose the pressure function p = p(p) is strictly increasing with respect 
to p > 0, the density p and velocity v of the Newtonian fluid are continuous and bounded 
in (x, t) G Q x [0, T*] with T* < oo ; and satisfy inf p(x, t) > and 



x.t 



Vp G C([0,T*},H s (n)), v G C^TtlH'+^O)) 

with integer s > 2. Then there exist positive numbers e = eo(T*) and K = K(T*) such 
that for e < e the first-order system M.3\) with initial data (p, pv, Tu, T 2e )\ t =o has a unique 
classical solution (p e , p t v t , T-f, t|) satisfying 

(p* - p,p*v%Tlr*) e C([0,T^H s (n)) 

and 



sup ||[(p e ,pV,7f,7£) - (p,/w,T le ,r 2e )](-,t)|| s < K(T*)e 2 
te[o,r«] 



(3.2) 



Here f2 = M 3 or the three-dimensional torus [0, l] 3 , and we are using the standard notation 
for Sobolev spaces, defined in [2l El El HH H21 H31 HU - 
An obvious corollary of this theorem is 

Corollary 3.2. If p and v possess the properties assumed in Theorem \3. 1\ globally in time, 
then the time interval where (p e , p e i> e , r^, r 2 ) exists goes to infinity as e tends to zero. 
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To see the existence claim in Theorem 13. 1[ we recall from the previous section that the 
first-order system (jl.3p is symmetrizable hyperbolic. Thus, the local-in-time existence 
theory of regular solutions to initial-value problems of symmetrizable hyperbolic systems 
applies j6]. Fix e > 0. According to the local-in-time existence theory, there is a time 
interval [0,T] such that (11.31) with initial data (p, pv, Tu, r 2e )| t= o has a unique solution 
(p e , p e u e , rf, r|) satisfying 

(p e - p, pV, r[, r*) G C([0, T] , H S (Q)). 

Note that the range G\ of (p, pv , Tu, r 2t )(x, t) satisfies G\ CC G for inf p(x,t) > 0. For 

x,t 

G 2 C G satisfying G\ CC G 2 , we define 

T e = sup{T > : (p e - p, pV, rl r 2 e ) G H S (Q), (p e , pV, rf, r 2 e )(x, t) G G 2 }. 

Namely, [0,T e ) is the maximal time interval of if s (f2)-existence. Note that T e = T e (G 2 ) 
may tend to as e goes to 0. 

In order to show that T e > T*, we exploit the convergence-stability lemma [T2l [2] and 
only need to prove the error estimate in (13.21) for t G [0, min{T*, T 6 }). In this time interval, 
both (p e , ,tI,t 2 ) and (p, pv,Ti e ,r 2e ) are well defined, regular enough and take values 
in the compact set G 2 . 

4. Error Estimate 

The purpose of this section is to derive the error estimate in (13.21) for t G [0, min-fT!,,, T e }). 
To do this, we need some classical calculus inequalities in Sobolev spaces [U E] . 

Lemma 4.1. (i). For s > 2, H s = H S (Q) is an algebra. Namely, if f,g G H s , then 
fg G H s and, for all multi-indices a with \a\ < s, 

\\d:(f9)\\<C s \\f\\ s \\g\\ s . 

Here C s is a generic constant depending only on s. 

(ii). For s > 3, let f G H s and g G if* -1 . Then for all multi-indices a with \a\ < s, 
the commutator [d®, f]g = d^(fg) — fd^g G L 2 (Q) and 

\mj]g\\ <C a ||V/|| a _i||^|| s _i. 

(Hi). Let f(u) be a smooth function of u. Then for all multi-indices a with \a\ > 1 we 
have 

\\d2f(u)\\ < CdttlooJHio,, 

where C(|m|oo) is a constant depending on the maximum norm \u\oo of function u = u(x). 
Now we derive the error estimate. With the notation in Section 2: 

W E = {p% pV) T , = (Yi, T 2 ) T , W e = (p e , p e V e ) T , Z e = (n e! T 2t f '. 

we set 

Ei = w e — w e and E 2 = z e — z e . 
From (12. 2p and (13. ip we deduce that 

dtEt + Ejifj^) - fi(«; )h; + \ Ej CjE 2Xj = 0, 

(4.1) 

d t E 2 + \ J2j(9j(™ e ) ~ 9j(w € )) x . = -±SE 2 - d t (r u , r 2e ) T . 
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Let a be a multi-index with \a\ < s. Differentiating the two sides of the last equations 
with and setting 



E la = d%E 1 and E 2a = d%E 2 , 

we obtain 

d t E la + £ . Aj(w*)E la:Cj + i E, C^2 OT , = / Q , 

9^2a + 7 E, B^E^ = -jzSE 2a - d t <9£(r le , r 2e ) T + \g a , 

where A,(u>) = f jw (w),Bj(w) = g jw (w), 

r = Ej -i ; !«•')/•„,,, - Ej - 

<f = Ei B^)E lax . - Ej 9?G&(«; e ) - <?> £ )W 



(4.2) 



Set = rj ww (w, z) and = r) zz (w, z). Note that D(ix;)Aj(iu) is symmetric and the 

identity f]2 .3f) holds. Multiplying the first equation in (14.21) with Ef a D(w e ) and the second 
with summing up the two and integrating the resultant equality over Q gives 

i J n [Ef a D{w*)E la + E 2 r a HE 2a ]dx 

= -f J n ElHSE 2a dx + 2 / n Ef a D(w^)f a dx - 2 /„ E^Hd t d^(r le , r 2e ) T dx 
+f /„ £ 2 T a # ^ + /„ + E i D (w e )Aj (w e )] E la dx 

+! Ej Jn El a d x p{w^)CjE 2a dx (4.3) 

< -^ll^ll 2 + 0||i?2 Q || 2 + Ce 4 + C||£ lQ || 2 

+c(||r II 2 + Ibll 2 ) + IW^) + Ej^j^M^iMUII^iair 

+C\T ij d Xj D{w^)C j U\Eiaf. 
Here c and C are both generic positive constants, and we have used that \\d^d t (ru, r 2t ) \\ < 



Newtonian limit of Maxwell fluid flows 7 

Next we analyze the terms in the last two lines. For ||/ a ||, we use Lemma [4.11 and the 
boundedess of \\w ex . \\ s to estimate as follows. 

Ilrll = II EjlM^Elax, ~ - fj(w e )) Xj }\\ 

< Ej \\Mw*)E lax , - d%{Aj{w*)w%, - A 3 {w e )w eX] )\\ 

< Ej \\Mw*)E laXj - d^{A 3 (w^)E 1Xj + (Aj(w<) - Aj(w £ ))w eXj ) \\ 

< Ej \\[Mw*),d a }E 1Xj \\ + £\ II^K^M - Mw e )) WeXj )\\ 

< C-^- llV^C^^IU-ill^^JU.! + C s ^. - ^(^ e )|U]|tz7 e;c J| s (4.4) 

< c.Ej \\VM™Ms-i\\Ei\\ s + c s Ej \\VM™ e ) - v^-c^oiU-ill^ilU 

+c s E 3 \\M^)-M^)\\s 

< Csimu + c.(i + \\EkW.) Ej KK) - M^)\U 

<C s (l + ||^i|| s )||^i|| s . 
Similarly, we have 

\\g a \\ <C s (i + 1^1)1^1. (4.5) 

In addition, we have 

E,<l »(«■<)(; <ce,\w Xj I 

\d t D(w<) + Ej d Xj D(w*)Aj(w*)\ < CZj Kl + C|i 
For \wl\ we use the equation and z t = 0(e) to get 

Kl< cEjKil + c^EjKl 



< c + ce 3 \e, Xj \ + c±E 3 \E2 Xj \ 

< c + cii^n. + cin^iu 

Substituting (T4~4"]) -f H~T|) into fT4~3]) we arrive at 

| / n [Ef a D(w*)E la + ElHE 2a ]dx 

< ~M E ^f + Ce± + C(l + ||^||2 + WmDWErWl 
Summing up this inequality over all a with \a\ < s gives 

i E a L [El a D{w*)E la + ElHE 2a ]dx 

< -M E *\\ 2 S + + CO- + W\l + ill^y I|£iil2 

< + + c(i + + ^11^112 + cWEkW* 

< -^ll^e + ^ 4 + ^(i + ll^ll?)ll^ill^ 



(4.6) 



(4.7) 
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Integrating the last inequality from to t G [0, min{T e , T*}) and using the positive defi- 
niteness of D{w e ) and H, we get 

ll^ilL 2 + liml <Ce' + C f\l + WE^mH = (4.8) 

Jo 

Obviously, 0(0) = Ce 4 and 

<j ) >(t) = C(l + \\E 1 \\ 2 s )\\E 1 \\ 2 s <C<j ) (l + <!>). 
Applying the nonlinear Gronwall-type inequality in [11] to the last inequality yields 

<f>(t) < e CT * 

for t G [0, min{T e , T*}) if we choose e so small that (f)(0) = Ce A < e~ CT * . Finally, we apply 
the standard Gronwall inequality to (14. 8 j) to obtain 

M + ||^2|L 2 <<P(t)<Ce A e CT *. 

This completes the proof. 
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